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ABSTRACT 

We report the first observations of the clustering environment of damped Lyman a absorption systems 
at z ~ 4. Color selection and photometric redshifts were used to select 44 candidate Lyman-break galaxies 
brighter than Iab ~ 25.5 from deep BRI images of the 35 arcmin^ field containing the quasar BR 
0951—04. Multislit spectroscopy of 35 candidate galaxies was performed and 8 of these candidates have 
been confirmed as z > 3.5 Lyman-break galaxies. With only BRI photometry, the photometric redshifts 
are quite accurate for the spectroscopically confirmed galaxies but have a high rate of misclassification 
due to color degeneracies between Lyman-break galaxies and low-redshift ellipticals. Both of the z > 3.5 
galaxies found within 15" of the quasar line-of-sight appear to be causing absorption systems in the 
quasar spectrum. We use a battery of statistical tests to look for clustering in the redshift histogram 
of the z > 3.5 galaxies but do not find measurable clustering of these Lyman-break galaxies with the 
damped Lyman a absorbers. With a larger sample of galaxies, our method should determine the cross- 
correlation between these objects, which probes the bias and hence the mass of the damped Lyman a 
absorbers. 

Subject headings: galaxies: high-rcdshift, cosmology: observational 



1. INTRODUCTION 

Galaxy formation is one of the great unsolved problems 
in astrophysics. We have clear evidence that many galax- 
ies were in place at redshifts of 3 — 4 but little knowledge 
of how or w hen they formed. Us ing Lyman-break imaging 
techniques, Steidel et al. (1996 ) discovered a large sam- 
ple of z ~ 3 galaxies with rest-frame ultraviolet luminosi- 
ties Ljjv > lOL* and comoving number density compa- 
rable to that of galaxies today. — These Lyman-break 
galaxies (LBGc) exhibit eurprieingl}^ ctrong cluetoring for 



evolution models in which the DLAs are more rarely oc- 
curring massive objects. Determining the mass of indi- 
vidual DLAs is also critical for testing models of galaxy 
formation within the CDM paradigm. This determination 
cannot be performed using the absorption characteristics 
of the DLAs, as the comoving density of DLAs cannot be 
determined without knowing their typical cross-sectional 
area. One way to distinguish among m odels for the DLAs 
is to raeasu r e the kinematics of the gas flProehaska & Wolfe 



1997, 1998; Wolfe fc Prochaska 2000a ,p:|). This paper de 



scribes an independent method which seeks to determine 
the mass of the DLAs by investigating the extent to which 
they cluster with LBGs. Since DLAs are predicted by hi- 
erarchical cosmologies to be rather weak in emission, an d 
observation is consistent with this ( Steidel et al. 1995 ), 
it may be easier to constrain their mass using large-scale 
stru cture techni ques than by searching for their emission 
(see [Wolfe 19"93| ). 

Section |^ describes the theoretical basis for our investi- 
gation of the cross-correlation of DLAs and LBGs. Imag- 
ing observations and photometric data reduction are de- 
scribed in and our multislit spectroscopy is presented 
in We discuss the search for emission from the DLAs 
themselves in ^||, search for evidence of clustering with 
nearby LBGs in §^ and discuss our results in §0. 

2. LARGE-SCALE STRUCTURE AT HIGH REDSHIFT 

Hierarchical models of structure formation predict that 
the most massive galaxies tend to form in regions of un- 
usually high density, whereas low-mass galaxies are more 
uniformly distributed in space. This leads to an enhance- 
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z ~ 3, with overdensities 5p/ p > 3 and comoving cor rela- 
tion length ro « 5h~^ Mpc ( [Adelberger et al. 1998| ). In 
hierarchical models of structure formation such as ACDM, 
the LBGs will undergo significant merging and should be 
even more clustered, rare massive galaxies today. 

The damped Lyman a absorption systems (DLAs), on 
the other hand, are th ought to be the progenitors of typica l 
present-day galaxies ( Wolfe et al. 1986 ; Kauffmann 1996|) . 
We therefore must learn about the mass and environment 
of high-redshift DLAs if we want to understand the forma- 
tion of typical galaxies like the Milky Way. The damped 
Lyman a absorbers are dense clouds of neutral hydrogen 
seen in absorption along the line of sight to distant quasars. 
They appear to be drawn from the bulk of the protogalac- 
tic mass distribution; this conclusion is supported by the 
agreement between the comoving mass density of neutral 
gas at z > 2 and t he mass density o f visible stars in cur- 
rent galaxy disks (Wolfe et al. 1995). However, the total 
gas content cannot distinguish CDM scenarios in which 
DLAs are numerous low-mass protogalaxies from passive 
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ment in the clustering of high-mass galaxies, referred to as 
bias, since they typically form near each other in rare over- 
dense regions. Clustering thus provides a way to probe the 
mass distribution of DLAs; if there is a significant over- 
density of LBGs near damped Lyman a absorbers it in- 
dicates that the DLAs reside in massive dark halos. If, 
on the other hand, the DLAs are preferentially found in 
regions with an underdensity of LBGs, it would argue that 
some environmental effect, perha ps ionizing radiation from 
massive stars in the LBGs (see Bteidel et al. 2000), pre- 
vents the occurrence of large column densities of neutral 
hydro gen in ovcrdensc regions. The expected result in hi- 



using SExtractor (Bertin & Arnouts 1996) and then a cus- 
tom code was used to extract photometry from the final 
B and I images using the exact apertures selected in R. 
Prochaska et al. (in preparation) describe our data reduc- 
tion methods in greater detail and make our photometric 
catalog and images public. 

We modified the G5Ri color selection criteria of S99 for 
our Johnson-Cousins BRI filter set.^ Objects located in 
the region of color-color space bounded by {B—R)ab > 
2.0, {B-R)ab > 1.2(R-1)ab + 1-6, and (R-I)ab < 1-0, 
have a high probability of being Lyman-break galaxie s at 



erarcl ical cosmologies is intermediate; DLAs are predicted 



to be low-mass objects and should be found in both the 
(proto)clusters and the field with typically only a small 
enhancement of LBGs nearby. 

While the cross-correlation of ga laxies and clust ers has 
been studied at low redshift (e.g. Mo et al. 1993), these 
objects can also be studied through their autocorrelation. 
The inevitable sparse sampling of DLAs and other absorp- 
tion systems caused by only being able to detect them on 
the line of sight towards quasars or bright galaxies makes 
it very difficult to measure their autocorrelation. Hence, 
the cross-correlation is the best, and perhaps only, way to 
study the bias of DLAs. Another advantage is that the 
cross- correlation of DLAs and LBGs can be established 
witho it detecting emission from the galaxy causing the 
damped Lyman a absorption. 



z > 3.5. We have also used photometric redshifts (see Ya, 
hata et al. 2000| ) and have found at least one LBG that 



would have been missed using only the two-color region of 
S99. Almost all of the color-selected candidates are identi- 
fied by the photometric redshift technique as having a high 
probability of being at 3.5 < z < 4.5. The photometric 
redshift method accounts for the individual photometric 
uncertainties in the B,R, and I fluxes of each object and 
compares the object fluxes with a set of template spectra 
rather than a single two-color region. Carefully designed 
seven-color photometric reds hift schemes yield a scatter of 
Az = 0.2 at these redshifts ( |Hogg et al. 1998| ). It should 
therefore be impossible for photometric redshift selection 
to create a redshift distribution as narrow as the Az ~ 0.04 



_(^2000 km/s) redshift spikes seen by Adelberger et al 
(199^ ) 



The probability of flnding an LBG in a small volume 
dV a distance r from a known DLA is enhanced versus 
the average density of LBGs, 



dPLBG = [1 + ^DLir)]pLBGdV 



(1) 



where S^dl is the cross-correlation function of DLAs and 
LBGs. On scales where linear bias is a good model, the 
cross-correlation of any two classes of objects is determined 
by multiplying the dark matter correlation function ^(r) 
by the product of the bias factors of each class: 



(2) 



The bias of a set of collapsed objects of a given mass can 
be d etermined using the extended Press-Schechter method 
(see Mo et al. 1996 J ). The autocorrelatio n of LBGs 
is known at z = 3 ( [Adelberger et al. 1998]) but h as not yet 
been determined aX z — A ( pteidel et al. 1999 , hereafter 



We visually inspected all candidate galaxies in our im- 
ages and rejected roughly half of the initial candidates 
for having suspect photometry. Many of the rejected ob- 
jects lied near the edge of our imaging field in regions of 
lower signal-to-noise. It is no surprise that such a high 
fraction of the initial candidates had suspect photome- 
try, as the color selection region is sparsely populated 
and therefore the small fraction of low-redshift galaxies 
scattered into the color-selection region by photometric er- 
rors can easily comprise a large fraction of the candidates. 
Bright (I< 23) objects with clearly non-stellar morphology 
were rejected as low-redshift interlopers, but bright stel- 
lar objects were retained in case they were new quasars or 
lensing-magnified LBGs; most LBGs will look similar to 
point sources at the resolution of our images. 

4. SPECTROSCOPY AND REDSHIFT DETERMINATION 

We performed longslit spectroscopy of two objects near 
BR 0951—04 and multislit spectroscopy of this field with 
four different slitmasks, all with 1.5" width slitlets and 
the 150 line/mm grating (30 A resolution) on Keck-LRIS. 
Table |l| lists details of our observations. Each slitmask 
targeted about 15 candidate objects. The entire LRIS 
imaging field was covered using 2 slitmasks (termed left 
and right) that overlap in the region near the quasar. 
These masks were redesigned for the March 2000 run 
to eliminate slits on objects whose redshifts were al- 
ready determined and to add slits on additional candidate 
LBGs. We used a GG495 dichroic filter to block second 
order light in December 1999 but decided against this 
in the subsequent run because wavelengths below 4950 
A can be used to check for the presence of the Lyman 
limit. The spectroscopic data were reduced using stan- 
dard IRAF techniques called by a customized version of 

®Our filter set corresponds closely to their custom set except for an overall 200A (A^ ~ 0.2) shift bluewards in B and R, allowing the 
approximation (B— R)yifl =(G — 5R). Because our I filter is close to theirs in effective wavelength, our (R— I)^^ colors should be roughly 1.6 
times their (SR — i) colors. 
'^Available at 



S99), probably due to the increased difficulty of obtaining 
redshifts at z ~ 4. The bias of other high-redshift objects 
such as quasars, radio galaxies, or Lyman a emitters can 
also be determined using the cross-correlation function. 

3. IMAGING AND PHOTOMETRY 

We took deep BRI imag es of the field around the z = 
4.37 quasar BR 0951-04 ( ^torrie-Lombardi et al. 1996D 
with the LRIS instrument (Oke et al. 1995) at the Keck 
Observatory. This field was chosen because there are two 
DLAs along the line o f sight towards BR 0951—04, a t 
z = 3.86 and z = 4.20 ( ptorrie-Lombardi fc Wolfe 2000] ). 
Table |l] lists the exposure time and final depth in each fil- 
ter. Images were reduced using a mix of standard IRAF 
routines and customized code. The final images cover a 
5' X 7' field. Photometry was performed on the R image 



http:/ /zwolfkinder .jpl.nasa.gov/~stern/homepage/bogus. html 
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the BOGUS code^. The functions response (halogen lamp 
spectrum fitting), baekground (sky subtraction) and apall 
(spectrum extraction) were performed interactively fol- 
lowing the multislit data reduction cookbook available at 



http: / /mamacass. ucsd.edu/people/gawiser/cookbooks. html 

Kedshifts were determined interactively using a cus- 
tomized PGPLOT routine which compared the measured 
spectra with template s of expected abso r ption (and emis- 
sion) lines for LBGs ( ^teidel et al. 2000 ; Lowenthal et al_ 
1997|]and E+A galaxies (|Zabludoff et al. 1996|). GamsCT 



et ah (in preparation) describe our spectroscopic data 
reduction in greater detail. 

We obtained spectra of 35 candidate galaxies in the field 
of BR 0951-04 (see Table ||). Eight candidates turn out to 
be galaxies at z > 3.5, thirteen appear to be low-redshift 
interlopers (typically ellipticals at z ~ 0.5) and fourteen 
have inconclusive spectra. Our color-selected sample with 
Iab < 25 is directly comparable to the objects observed 
by S99. We find 0.5 candidates per sq. arcmin, twice the 
candidate density of S99; the difference is unlikely to be 
caused solely by Poisson fluctuationsgiven that there are 
19 objects in this category in Table ||. Our rate of deter- 
mining redshifts for objects observed spectroscopically is 
comparable to S99, and the fraction of candidates iden- 
tified as interlopers, i.e. at z < 3.5, is higher but very 
sensitive to small- number statistics at present. We will re- 
visit this comparison once a larger sample of BRI-selected 
galaxies is available. 

Figure ll| shows the two-color diagram of all 2254 ob- 
jects with < 25.5 detected in the field with the status 
of LBG candidates shown. Table || lists photometry and 
photometric redshift information on each z > 3.5 galaxy 
with its best-fit redshift and uncertainty. The photometric 
redshifts based on BRI alone are impressively accurate for 
the confirmed high-redshift galaxies, with a dispersion of 
a = 0.18, including three cases where Zphot differs from the 
spectroscopic redshift by only 0.01. However, adding in the 
catastrophic failures where z ~ 0.5 galaxies were assigned 
photometric redshifts Zphot — 4 would drastically increase 
this dispersion. Adding near-infrared (JHK) or optical Z- 
band photometry would be quite helpful in discriminating 
against these intrinsically-red low-redshift interlopers. 

The increased fraction of interlopers among candidates 
selected primarily by photometric redshift confirms that 
the S99 color selection region is the most efficient place 
to find LBGs. However, we differ from S99 in one major 
respect. We do not assume that the spectra with insuffi- 
cient signal-to-noise for redshift determination are drawn 
from the same population of objects as those with iden- 
tified redshifts. Thus there is a factor of two uncertainty 
in the number abundance of color-selected z > 3.5 galax- 
ies, both here and in S99, caused by almost half of the 
candidates lacking redshift identifications. Extrapolating 
from the observed color distribution of LBGs at z = 3, 
S99 predict that roughly half of the z = 4 LBGs will lie 
outside of the color selection region. The failure of photo- 
metric redshifts to find these LBGs implies either that this 
incompleteness has been overestimated or that the photo- 
metric redshift technique needs to be modified to better 
identify LBGs with unusual colors. 

Figure ^ shows the rest-frame spectra of the confirmed 
z > 3.5 galaxies. These spectra show the typical mix of 
Lyman a emission (G4,G6,G7) and Lyman a absorption 
(G1,G2,G8) as well as two vague suggestions of emission 
within a broad absorption trough (G3,G5). The signal- 



to-noise is generally quite poor, illustrating why 1=25 is 
usually considered the limit for ground-based optical spec- 
troscopy. 

5. SEARCH FOR EMISSION FROM DAMPED LYMAN a 
ABSORBERS 

At z ~ 4, Lyman-break galaxies found within a 10" ra- 
dius of the quasar may be the source of damped Lyman a 
absorption. This radius corresponds to the maximum im- 
pact parameter (~ 100 kpc) inferred fro m absorption line 
statistics in the passive evolution mo del ( ^torrie-Lombardi 



fc Wolfe 2000| ; [Bunker et al. 1999] ). I f DLAs repr esent 



formed disk galaxies as postulated by Wolfe (1997), we 
should be a ble to detect emissi on from up to 20% of them 
at R< 25 (Mailer et al. 2000). However, in hierarchical 
models of structure formation, a much higher fraction of 
DLAs are expected to be too dim a nd/or too close to the 



quasa r to yield detectable emission (Haehnelt et al. 1998 



2000). Our survey hopes to differentiate between these 



hypotheses; if emission from DLAs is found repeatedly at 
large impact parameters or high luminosity it will repre- 
sent a serious conflict with the predictions of hierarchical 
cosmology. 

We performed a complete survey for possible emission 
from DLAs on the line of sight toward the quasar down 
to Iab — 26 at angular distances up to 10", although the 
PSF of the quasar hampers our survey at angles less than 
2" where hierarchical cosmology predicts that the emis- 
sion from most DLAs will be found. We identified three 
objects within 10" of the QSO line-of-sight that had a 
reasonable chance of being at high redshift (our selection 
criteria were less stringent in this region). One of these, 
Gl, meets the color selection criteria but has Iab = 25.7 
and appears to have a redshift of 3.69. Although Gl is 10" 
from the QSO, it may be responsible for the z = 3.703 ab- 
sorption sy stem seen in SilV and CIV ( Btorrie-Lombardi 
et al. 1996| ). The second nearby object appears to be a 



low-redshift interloper, and the third has inconclusive red- 
shift. Another color-selected galaxy, G2 at z = 3.81, lies 
15" (~ 200 kpc) from the QSO and may be responsible 
for the z = 3.818 absorption system seen in Lyman a 
and CIV in our Keck-HIR ES spectrum of BR 0951-04 
( Prochaska fc Wolfe 1999| ). It will be interesting to de- 



termine the fraction of galaxies detected in both emission 
and absorption as a function of distance from the quasar 
line-of-sight. We have not detected emission physically 
associated with either of the damped Lyman a absorbers 
towards BR 0951—04 spectroscopically, but all of the mod- 
els predict that a larger sample is needed before we should 
expect to find an LBG physically associated with a DLA. 

6. CLUSTERING ANALYSIS 

Figure ^ shows the redshift histogram of objects in the 
field ofBR 0951-04, with a binsize of Az = 0.05 (~ 10/1"^ 
Mpc) chosen because most of our redshift uncertainties are 
smaller than this value. We chose to place the z = 3.86 
DLA in the center of a bin in order to maximize the chance 
of finding an overdensity (or underdensity) of LBGs near 
the DLA if one exists. However, the results described be- 
low are not altered by shifting the bin locations. 

There are eight confirmed LBGs in the BR 0951-04 
field, so the search for clustering is plagued by small num- 
ber statistics. Before trying to determine the amplitude of 
a cross-correlation function, it makes sense to ask if any 
clustering has been observed that is inconsistent with pois- 
son fluctuations about the expected redshift distribution 
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of LBGs. The expected redshift distribution refers to the 
multipUcation of the true redshift distribution of LBGs on 
an average hne of sight times the observational selection 
function. While the observational selection function can 
be simulated (see S99), the true redshift distribution of 
LBGs is unknown, so it is necessary to determine the ex- 
pected redshift distribution empirically. We have not yet 
observed a sufficient number of LBGs to determine how 
our expected redshift distribution differs from that of S99, 
who found < z >= 4.13, cr = 0.26. We do expect a shift 
of 0.2 in redshift, so we choose to bracket the possibilities 
by showing three possible expected redshift distributions 
in Figure y. The mean of all three is z = 3.9; the dashed 
one is a gaussian with a = 0.3 and is considered the most 
likely, the dot-dashed one is a gaussian with a = 0.2 and 
is considered the narrowest possible, and the flat dotted 
one is the broadest possible (a uniform distribution from 
z = 3.3 to z = 4.5). The presence of a bin in our redshift 
histogram that contains three LBGs is only a weak indica- 
tion of autocorrelation of the LBGs, because Poisson fluc- 
tuations of any of the three expected redshift distributions 
generate an appearance of clustering at least this strong 
over 10% of the time in our Monte Carlo simulations.^ A 
Kolmogorov-Smirnov test of the unbinned cumulative red- 
shift distribution versus the three expected distributions 
fails to find any inconsistencies at the 95% confidence level. 
A similar analysis of the reported redshifts of S99 indicates 
that the redshift histograms of all 7 of their fields are con- 
sistent with Poisson fluctuations of their expected redshift 
distribution. 

We are interested primarily in the clustering environ- 
ment of the DLAs. This allows us to select the bins known 
a priori to contain a DLA and to ask if a statistically signif- 
icant overdensity (or deficit) of LBGs has been found there. 
Finding zero LBGs in the bin containing the z = 3.86 DLA 
should occur at least 40% of the time for any of the ex- 
pected redshift distributions. Finding zero LBGs in the 
bins containing the z — 4.20 DLA or z = 4.37 quasar is 
at least as common. Even the apparent triple coincidence 
of finding no LBGs in any of those bins occurs in at least 
30% of our Monte Carlo simulations, so we do not have 
good evidence for an anticorrelation of DLAs and LBGs. 
Our results provide an upper limit on the amplitude of the 
DLA-LBG cross-correlation function, but it is not a partic- 
ularly interesting upper limit, as even the observed strong 
autocorrelation function of z = 3 LBGs would be diffi- 
cult to detect with so few objects. A vast improvement in 
statistics will be required to measure the cross-correlation 
of LBGs and DLAs. 

The minimum redshift uncertainty of our LBGs is Az — 
0.01, so correlations on comoving scales of a few Mpc and 
smaller can only be probed via angular clustering. While 
8 LBGs and 2 DLAs are insufficient to measure the angu- 
lar cross-correlation precisely, we can again test whether 
the observed distribution is consistent with being random. 
The presence of 1 out of 7 LBGs from our Iab < 25.5 
sample within 15" of the QSO line-of-sight is inconsistent 
with a random distribution at the 95% confidence level. ^ 
Since G2 is not particularly close to the QSO or DLAs 



in redshift this angular overdensity appears to be a coin- 
cidence rather than the projection of a three-dimensional 
structure. 

7. DISCUSSION 

We have not found significant evidence for autocorrela- 
tion amongst the 8 LBGs at z > 3.5 in our current sample. 
It is therefore not surprising that we find no evidence for 
a cross-correlation between LBGs and damped Lyman a 
absorbers, as the latter should be more difficult to detect 
if the DLAs are indeed less massive than the LBGs. Our 
results at z ~ 4 are comparable to the limited studies of 
the clustering environment of dampe d Lyman a systems 



that have been performed at z ~ 3. Steidel et al. (1996) 
found only one out of eleven Lyman-break galaxies in the 
field of QOOOO-263 to be within Az = 0.04 of the z = 3.39 
DLA. [Ellison et al. (2000D found two out of 27 LBGS in the 
field containing Q0201+1120to be within Az = 0.02 of the 
z — 3.39 DLA. This failure to detect clustering of LBGs 
with DLAs at z ~ 3 stands in contrast to the detection 



of strong a utocorrelation of LBGs at z ~ 3 (Adelberge: 
et al. 199^ ), although the latter was determined from a 
much larger sample of galaxies. 

In order to obtain better statistics in the future, we have 
begun a systematic survey for Lyman-break galaxies asso- 
ciated with damped Lyman a absorbers at z ~ 3. The 
increased sky density of LBGs bright enough for spec- 
troscopic identification at z ~ 3 and reduced sky noise 
near 5000 A should allow us to obtain significantly im- 
proved statistics. The coincidence of LBG overdensity and 
damped Lyman a redshifts in a number of fields would ar- 
gue in favor of strong biasing and hence large mass of the 
damped Lyman a absorbers, which is counter to the pre- 
dictions of the ACDM model. 

We see no positive indication of clustering of LBGs with 
the quasar BR 0951—04, but this quasar is at too high a 
redshift for us to be sensitive to an LBG-quasar cross- 
correlation except in the unlikely case that our expected 
redshift distribution is uniform. Observing fields with a 
quasar at lower redshift will allow us to measure the cross- 
correlation of quasars and LBGs, which will test the belief 
that quasars are highly biased markers of protoclusters 
(Djorgovski 1997). With better statistics, this observa- 
tional program will measure the bias of DLAs, quasars, 
and Lyman-break galaxies, revealing much about the typ- 
ical mass and environment of these high-redshift objects. 
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"Because this bin centered at z = 3.71 is below the mean redshift, the result is more significant for the gaussian expected redshift distribu- 
tions, but this could easily be weakened if the true mean of the expected redshift distribution is 2 = 3.8. We will revisit this issue when the 
expected redshift distribution is better determined by observations in additional fields. 

^Gl was selected only because it lies within 10" of the quasar and was therefore excluded from this analysis. 
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Table 1 

Observing log for imaging and spectroscopy in the field of BR 0951—04. 



Date 


Observing Mode 


Exposure Time 


Seeing 


la sky 


3(7 source detection 






(s) 


(") 


AB" mag/sq. " 


AB mag/FWHM^ 


1999 January 13 


B imaging 


4500 


0.85 


29.7 


28.7 




R imaging 


1620 


0.75 


29.0 


28.1 




I imaging 


1560 


0.63 


28.1 


27.4 


1999 March 22 


longsht 


7200 


1.2 






1999 December 6 


multisht-left 


7200 


1.0 








multisht-right 


3600 


1.3 






2000 March 9 


multisht-left 


3600 


1.1 








multisht-right 


7200 


1.2 







" An AB magnitude of zero corresponds to 3.64 x 
converted to the Johnson-Cousins system using Bjc 



10-23 -y^ ^-2 AB95 magnitudes ([Fukugita et al. 1996| ) can be 

= Bab + 0.15, Rjc = Rab - 0.18, and Ijc = Iab - 0.43. 
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Table 2 

Spectroscopic results for candidate z > 3.5 galaxies. 



Category: 


I 


II 


III 


IV 


V 


all 


z > 3.5 


3 


1 


3 





1 


8 


z < 3.5 


3 


6 


3 


1 





13 


uncertain 


8 


1 


4 


1 





14 


not observed 


5 





5 








10 


Total 


19 


8 


15 


2 


1 


45 



I: meets two-color criteria, Iab <25 

II: selected based on photometric redshift, Iab <25 

III: meets two-color criteria, 25 <Iab < 25.5 

IV: selected based on photometric redshift, 25 <Iab < 25.5 

V: within 10" of quasar, 25.5<Iab <26.0 



Table 3 

SpECTROSCOPICALLY confirmed Z > 3.5 GALAXIES. 



Galaxy: 


z 




^AB 


(B-R)^B 


{R-1)ab 


7 y.^ 

^proo 


rv C 

^phot 


Gl 


3.69 


-0.05 


25.7 


2.2 


0.0 


0.74 


3.78 


G2 


3.81 


0.01 


24.6 


2.8 


0.0 


1.00 


3.82 


G3 


3.72 


+0.03 


24.2 


2.8 


0.0 


1.00 


3.85 


G4 


3.57 


0.01 


24.8 


1.6^^ 


-0.2 


0.82 


3.62 


G5 


3.75 


-h0.05 


25.3 


2.6 


-0.1 


0.91 


3.76 


G6 


4.41 


0.01 


24.9 


3.2 


0.0 


0.99 


3.93 


G7 


3.72 


0.01 


25.2 


2.3 


-0.2 


1.00 


3.71 


G8 


4.14 


+0.09 


25.5 


3.6 


0.2 


0.90 


4.08 



° A plus or minus sign indicates that a specific other redshift, z +(t^ is also a good fit. Otherwise, the errors are approxi- 
mately symmetrical and roughly correspond to 95% confidence. 

^ Defined as the integrated probability that the redshift is in the range 3.5 < z < 4.5. The method corresponds to 
assuming a uniform prior as a function of redshift but then optimizing the likelihood at each redshift over a set of possible 
template spectra. 

The redshift which maximizes the photometric redshift likelihood function. 

G4 falls outside of the color selection region and was included on the basis of its photometric redshift. The very blue 
R — I color makes it fall directly on the evolutionary track of starburst galaxies at 2: ~ 3.6 which are at too low redshift 
to drop out of the B filter as fully as the color criteria require. 
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Fig. 1. — Two-color diagram of objects imaged in field of BR 0951—04 with \ab < 25.5 (with the addition of object Gl). Objects 
that were selected as candidate LBGs based on their colors and/or photometric redshifts have been assigned large symbols as 
identified in the upper left corner. Objects that fall in the color selection region but are represented only by small points were 
found by visual inspection to have suspect photometry and were therefore not considered as candidate LBGs. 
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Fig. 2. — Spectra of 8 confirmed z > 3.5 galaxies in the field of BR 0951—04 after being smoothed by the resolution element. 
The dotted spectrum shows the \g uncertainty. Spectra are plotted versus inferred rest-frame wavelength, with dotted vertical 
lines showing expected absorption/emission features. Flux is shown in units of /xJy (1/iJy = 10~ erg cm" -1 s"^ Hz-i). Dashed 
vertical lines show sky lines (at 5577, 5892, 6300, 6364, 6840, 6940, 7244, 7276, 7340, and 7475 A ) that are too strong to be 
subtracted well and can be used as a rough guide to the observed wavelengths. Observed wavelengths from 4000 to 7500 A are 
shown; the signal-to-noise degrades considerably below that due to reduced grating throughput and CCD quantum efficiency and 
above that due to sky emission bands. Object G6 was observed with a GG495 dichroic, so only wavelengths above 4900 A are 
shown. 



Damped Ly a Absorbers at z ~ 4 
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Fig. 3. — Histogram of redshifts of LBGs in field of BR 0951—04. Redshifts of DLAs and quasar are shown by vertical dotted 
lines. The range of possible expected redshift distributions described in the text is also plotted; the dot-dashed line has mean 
z = 3.9, a — 0.2, the dashed line has mean z = 3.9, a — 0.3, and the dotted line is a uniform distribution from z = 3.3 to z = 4.5. 
All are normalized to produce 8 galaxies. 



